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Abstract 

In this paper we study the detectability of 7-rays from dark matter annihilation in the subhalos 
of the Milky Way by the satellite-based experiments, EGRET and GLAST. We work in the frame 
of supersymmetric extension of the standard model and assume the lightest neutralino being the 
dark matter particles. Based on the N-body simulation of the evolution of dark matter subhalos 
we first calculate the average intensity distribution of this new class of 7-ray sources by neutralino 
annihilation. It is possible to detect these 7-ray sources by EGRET and GLAST. Conversely, if 
these sources are not detected the nature of the dark matter particls will be constrained by these 
experiments, depending, however, on the uncertainties of the subhalo profile. 
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I. INTRODUCTION 

In the third EGRET catalog |l| unidentified 7-ray sources accounted for more than a half 
of the discrete sources detected by EGRET. Despite a great effort, most of them can not be 
associated with the known sources detected at low energies up to now. Most of the efforts 
in identifying these sources focused on the Galactic counterparts, such as young pulsars 
j^], microquasars j^], supernova remnants At the same time, multiwavelength searches 
continue to look for counterparts in the 7-ray sources, such as blazars, supernova remnants 
and pulsars [5J. There are recent theoretical works trying to explain the unidentified sources 
as annihilating dark matter clumps within the Milky Way (MW) {(]]. These efforts try to 
solve the nature of dark matter. 

The existence of cosmological dark matter has been firmly established by a multitude of 
observations, such as the observations of the rotation curves in spiral galaxies and velocity 
dispersion in elliptical galaxies, the X-ray emission and peculiar velocities of galaxies in 
the clusters of galaxies, the weak lensing effects, all indicating much steeper gravitational 
potentials than those inferred from the luminous matter. However, the nature of the non- 
baryonic dark matter is still unknown and remains one of the most outstanding puzzles in 
particle physics and cosmology. 

Among a large amount of theoretical candidates, the most attractive scenario involves 
the weakly interacting massive particles (WIMPs). An appealing idea is that the WIMPs 
form the thermal relics of the early universe and naturally give rise to the relic abundance in 
the range of the observed values. The WIMPs are well motived theoretically by the physics 
beyond the standard model to solve the hierarchical problem between the weak scale and the 
Planck scale. In particular, the minimal supersymmetric extension of the standard model 
(MSSM) provides an excellent WIMP candidate as the lightest supersymmetric particle, 
usually the lightest neutralino, which are stable due to R-parity conservation ?]]. The cos- 
mological constraints on the supersymmetric (SUSY) parameter space have been extensively 
studied in the literature |8|]. 

The WIMPS can be detected on the present running or future experiments, either directly 
by measuring the recoil energy when WIMP scatters off the detector nuclei [9| or indirectly 
by observing the annihilation products of the WIMPs, such as the antiprotons, positrons, 
7-rays or neutrinos [10|] . The rate of the WIMP annihilation is proportional to the number 
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density square of the dark matter particles. Therefore the searches for the annihilation 
signals should aim at the re gion s with high matter densities, such as at the galactic center 



ll| or the nearby subhalos 
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14j . The existence of a wealth of subhalos throughout 



the galaxy halos is a generic prediction of the CDM paradigm of structure formation in the 
Universe. High resolution simulations show that for CDM scenario the large scale structure 
forms hierarchically by continuous merging of smaller halos and as the remnants of the 
merging process about 10% to 15% of the total mass of the halo is in the form of subhalos 
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22] . At the center of the subhalos there are high mass densities 



and therefore they provide good sites for the search of WIMP annihilation products. 



However, the analysis in Refs. jf], [23j shows that it seems most of the EGRET uniden- 
tified 7-ray sources are not produced by the dark matter annihilation. In jg] by comparing 
the cumulative luminosity function of subhalos and EGRET unidentified sources the au- 
thors found at most 26 ± 11 unidentified sources are possibly subhalos. However, location 
coincidence and variability cut even exclude further of these candidates jf|. In 23] assuming 
a similar population of subhalos between the Milky Way and M31 and using the upper limit 
of 7-rays from M31 the authors found it is highly unlikely that a large fraction of these 
unidentified EGRET sources can be from subhalos. Improving the upper limit may finally 
exclude this possibility [231 ] . 



The strong constraints on the possible detectable subhalos from these works [6|, [23| will 
constrain the properties of dark matter particles. In the present work we assume that 
neutralino forms dark matter and work in the frame of supersymmetric extension of the 
standard model. Assuming that none of the unidentified EGRET sources is from dark 
matter annihilation we will study how the properties of neutralino are constrained resorting 
to the numerical simulation result of dark matter clumps distribution. The non-detection 
of subhalos finally shows a constraint on the SUSY parameter space. 

The next generation satellite based experiment, GLAST 241 ] . will greatly improve the 



sensitivity of EGRET. More 7-ray sources will be detected by GLAST. Detectability o 



rays from the subhalos by GLAST has been studied in literature 13 
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3l|. Once such sources are detected, the follow-up study of the sources may measure the 



spectrum of the annihilated 7-rays or even detect the line emission 31] . These measurement 
will finally give strong implications on the properties of dark matter. However, there are also 
possibilities of null result for such searches, similar to the EGRET result or the null result 
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of direct detection even though the direct detection has continuously improved sensitivities. 
In this case we would like to study how supersymmetry would be constrained. 

It is well known that the flux of gamma rays from the neutralino annihilation in a clump 
is given by 

(av)dN f^o 2 1 (av)dN 1 f f , 2 , , 

where (erf) is annihilation cross section times relative velocity, 4^ is the differential flux in 
a single annihilation, m is the mass of the dark matter particle, d is the distance from the 
detector to the source. The flux depends on both the distribution of the dark matter p(r) 
and the particle nature of dark matter. 

In the next section we first give the intensity distribution of the subhalo 7-ray sources 
according to the N-body simulation results. Then in Sec. Ill we present the constraints on 
the SUSY parameter space from non-detection of the subhalo 7-ray sources by EGRET and 
GLAST. We finally give a conclusion in Sec. IV. 



II. DISTRIBUTION OF SUBHALOS AND THEIR 7-RAY INTENSITIES 



In this section we first present some simulation results about the subhalos and then we 
calculate the intensities of the MW subhalos as 7-ray sources by neutralino annihilation. 

N-body simulations show that the radial distribution of substructures is generally shal- 
lower than the density profile of the smooth background due to the tidal disruption of 
substructures which is most effective near the galactic center 32]. The relative number 



density of subhalos is approximately given by an isothermal profile with a core 



n(r) = 2n H {l + {r/r H ) 2 )- 1 



32] 



(2) 



where n.# is the relative number density at the scale radius r#, with r# being about 0.14 
times the halo virial radius r# = 0.14r V j r . The result given above agrees well with that in 



another recent simulation by Gao et al. 33]. 

The differential mass function of substructures has an approximate power law distribu- 



m 



with a 



-1.7 



32 



34). In Ref. |32|, |34| both the cluster 



tion, dn/dm 

and galaxy substructure cumulative mass functions are found to be an m -1 power law, 
n su b(m sub > m) oc m _1 , with no dependence on the mass of the parent halo. A slight 
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difference is found in a simulation by Gao et al. [33] that the cluster substructure is more 
abundant than galaxy substructure since the cluster forms later and more substructures 
have survived the tidal disruption. The mass function for both scales are well fitted by 
dn/dm oc m~ L9 . Taking the power index of the differential mass function greater than —2 
makes the fraction of the total mass enclosed in subhalos insensitive to the mass of the min- 



imal subhalo we take. The mass 
between 5 percent to 20 percent 



raction_of subhalos estimated in the literature is around 
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351 ] . In this work we will take the differential 



index of —1.9 and the mass fraction of substructures as 10 percent. 

To calculate the 7-ray intensities from the dark matter clumps we first realize a MW-like 
halo with a population of subhalos due to the subhalo distribution function from simulation, 
which is given above. The mass of the substructures are taken randomly between M min = 
1O 6 M , which is the lowest substructure mass the present simulations can resolve [3^], and 
the maximal mass M max . The maximal mass of substructures is taken to be 0.01M vir since 
the MW halo does not show recent mergers of satellites with masses larger then ~ 2x 1O 1O M . 
The 7-ray flux is quite insensitive to the minimum subhalo mass since the flux from a single 



subhalo scales with its mass 
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37]. 



However, due to the finite spatial resolution of the N-body simulations the distribution 
in Eq. ([2]) is an extrapolation of the subhalo distribution at large radius. The formula 
underestimates of the tidal effect which destroys most substructures near the Galactic Center 
(GC). We take the tidal effects into account under the "tidal approximation" , which assumes 
that all mass beyond the tidal radius is lost in a single orbit while keep its density profile 
inside the tidal radius intact. 

The tidal radius of the substructure is defined as the radius at which the tidal forces of 
the host exceeds the self gravity of the substructure. Assuming that both the host and the 
substructure gravitational potential are given by point masses and considering the centr ifug al 



force experienced by the substructure the tidal radius at the Jacobi limit is given by 



where r c is the distance of the substructure to the GC, M(r < r c ) refers to the mass within 
r c - 

The substructures with r tid < r s will be disrupted. The mass of a substructure is recalcu- 
lated by subtracting the mass beyond the tidal radius in realizing the MW-like halos. After 



taking the tidal effects into account we find the substructures near the GC are disrupted 
completely. 



A. concentration parameter 



We adopt both the NFW [39|, |40j and Moore |4ll . |42| profiles for the subhalos in our 



calculation, which can be written in a general form as 

9 (r/r s )~f[l + (r/r s ) a }(P-^)/ a, {> 

where p s and r s are the scale density and scale radius respectively (a, (3, 7) = (1,3, 1) and 
(1.5, 3, 1.5) are for the NFW and Moore profiles respectively. 

The free parameters p s and r s can be determined by the mass and concentration parameter 
of the subhalo. The concentration parameter is defined as 

c= , 5 

where r vir is the virial radius of the halo and r_ 2 is the radius at which the effective loga- 
rithmic slope of the profile is —2, i.e., J^(?" 2 p( r )) |r=r_ 2 = 0. For the NFW profile we have 
r s = r_ 2 , while for the Moore profile we have r s = r_ 2 /0.63. 

The concentration parameter is a crucial parameter in determining the 7 ray fluxes from 
subhalos. From the definition of the concentration parameter in Eq. ([5]) and the annihilation 
flux in Eq. (CQ) we can easily get that the annihilation flux from a clump is proportional 
to ~ Ap 2 s r z s ~ Arrive 3 with A a flat function of c. Therefore the annihilation flux is 
very sensitive to the concentration parameter. The concentration parameter is obtained by 
N-body simulations. However, due to the finite resolution of N-body simulations, numerical 
convergence has not been established, especially for the evolution of subhalos. Adopting 
different models the predicted detectable number of subhalos at GLAST (at 5cr for 1 year 
exposure) spans fan < 1 by Koushiappas et al. Q to at m0 st abont 300 by Baltz et 
al. [311 ] . Recently Pieri et al. tried to classify different cases by modeling the subhalos 
concentration parameter and found the detectable number of subhalos at GLAST ranges 



from ~ to ~ 40 



271 ] for different models they adopted. Another way for this kind of st udy 



is directly based on the simulation result, such as in the recent work by Diemand et al. 34] 



In this work we will adopt different analytic models or the fit formulas based on simulation 
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FIG. 1: Concentration parameter as a function of the virial mass calculated according to the 
Bullock model 43J]. The model parameters are taken as F = 0.015 and K = 4.4. The cosmology 
parameters are taken as CIm = 0.3, $7 a = 0.7, J7^/i 2 = 0.02, h = 0.7, <7g = 0.9 with three 
generations of massless neutrinos and a standard scale invariant primordial spectrum. Both the 
median and the ±lcr values of the concentration parameters are plotted. 

results about subhalos in the literature to discuss their detectability. Uncertainties of the 
simulation results are thus included in our study. 

We first introduce a semi-analytic model to determine concentration parameter given by 
Bullock et al. [43j which is built based on their simulation result. In the model, at an 
epoch of redshift z c a typical collapsing mass M*{z c ) is defined by cr[M*(z)] = S sc (z), where 
cr[M*(,2)] is the linear rms density fluctuation on the comoving scale encompassing a mass 
M*, 5 SC is the critical overdensity for collapsing at the spherical collapse model. The model 
assumes the typical collapsing mass is related to a fixed fraction of the virial mass of a halo 
M*(z c ) = FM vir . The concentration parameter of a halo with virial mass M vir at redshift z is 
then determined as c v i r (M v i r , z) = K^^. Both F and K are constants to fit the numerical 
simulations. A smaller M v i r corresponds to a smaller collapsing mass and early collapsing 
epoch when the Universe is denser and therefore a larger concentration parameter. Fig. 
CD plots the concentration parameter at z = as a function of the virial mass of a halo 

n 

according to the Bullock model [43J. 

From Fig. [T]we can see that between the masses 1O 6 M ~ 1O 1O M an experiential formula 
cvir oc reflect the simulation result accurately. We expect that this power law relation 
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FIG. 2: The cumulative number of subhalos as function of the integrated 7-ray fluxes n(> $ 7 ) for 
the Moore profile (left panel) and the NFW profile (right panel) within the solid angle of it. Sub- 
halos are taken as point sources within the angular resolution of 1°. The curves are corresponding 
to different simulation results, where 'subhalo' denotes the model for subhalos within a smooth 

n n 

halo [43]; 'Reed et al.' refers to the simulation results given by Reed et al in 14511 : 'Bullock et al.' 
uses the median c v i r — M„j r relation for distinct halos of the Bullock model in [43J]; 'ENS' refers to 



the result of Eke et al. 



44]. 



should be very well followed, since subhalos form early at the epoch when the Universe is 



dominated by matter with approximate power-law power spectrum of fluctuations 



43]. 



In the literature another widely adopted semi-analytical model for the concentration 



parameter is given by Eke, Navarro and Steinmetz (ENS) 44J. We also adopt the ENS model 



for the ACDM mode 
results by Reed et al 



with as = 0.9. The other two models we adopted are the simu 



ation 



451 ] and that by Bullock for the subhalos in dense environment [43J . We 
will show that these models predict very different annihilation fluxes. Especially for subhalos 
within the dense environment, simulation indicates it may have greater concentration than 
these of isolated halos 43j and therefore lead to larger annihilation flux. 



B. 7-ray intensity of the subhalos 

Once the profile parameters of each subhalo are determined and, in each realization of the 
MW-like halos, the distribution of subhalos is also known we can calculate the 7-ray fluxes 
from these subhalos. Then we can plot the number of sources as function of their intensities. 
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By realizing one hundred such MW sized halos we give the averaged number of sources as a 
function of their intensities. We have fixed the property of neutralino by requiring the 7-ray 
flux of 3.7 x 10~ 9 ph cm~ 2 s~ l sr~ l from the GC assuming a NFW profile. Fixing flux from 
the GC actually gives the relative intensities between the GC and each subhalo. 

Fig. Ogives the cumulative number of subhalos emitting 7-rays with fluxes greater than 
a value $ 7 . In the left panel we plot the results for the Moore profile while the right panel 
is for the NFW profile. From Fig. [2] we can easily read the expected number of detectable 
subhalos if the sensitivity of a detector is given with the same field of view. For example, 



the sensitivity of GLAST at 5a for 1 year exposure reaches 2 x 1Q~ 10 cm~ 2 s^ 1 sr~ l 46|. For 
comparison we adopt the same particle factor as Pieri et al. [2^], which is about 3 orders 
of magnitude than our particle factor adopted here and get the detectable subhalo number 
in NFW profile is from ~ to ~ 400. The maximal detectable number subhalos predicted 



by Pieri et al 27| is only about 40, however, which is given in a different scenario from the 



'subhalo' case here. Our result is consistent with the result by Baltz et al. [3JJ. Even for this 
scenario and with the optimistic particle factor we have only about 2 detectable subhalos at 
EGRET, whose 5a sensitivity can only reach ~ 10~ 8 cm -2 s -1 . 

Fig. [2] shows that there is a large discrepancy for predictions based on different models. 
Especially the 'subhalo' model gives much greater prediction. The reason is directly related 
with the large concentration parameters for the subhalo scenario. It should be noted that 
the other three models actually describe distinct halos with small masses. A qualitative 
simulation result about the concentration parameter is that it is determined by the halo 
collapse time, as shown in the Bullock model. The reason of large concentration parameter 
for 'subhalo' is that in dense environment halos tend to collapse much earlier [43]. Tidal 
stripping may also lead to stronger mass dependence of concentration on the subhalo mass. 
Another reason is that a high density environment likely leads to extreme collapse histo- 
ries of subhalos and frequent merger events which affect the final concentration (For more 



discussions see [431] ) 



C. detector sensitivities 

Before going on to the next section we first discuss the sensitivities of EGRET and 
GLAST. The detectability of a signal is defined by the ratio of the signal events to the fluc- 
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tuation of the background. Since the background follows the Poisson statistics, its fluctuation 
has the amplitude proportional to y/Ns- The significance of the detection is quantified by 

The signal events are given by 



n 7 



e A n [ A eff (E,6)(t>(E)dEdndT , (6) 

JE th ,An 



where = 0.68 is the fraction of signal events within the angular resolution of the in- 
strument and the integration is for the energies above the threshold energy E t h, within the 
angular resolution of the instrument AQ and for the observational time. Generally the ef- 
fective area A e ff is a function of energy and zenith angle. The <f)(E) is the flux of 7-rays 
from DM annihilation. We take E th = 1 GeV for both EGRET and GLAST. The EGRET 
has angular resolution of ~ 1° while GLAST has much better angular resolution of ~ 0.1°. 

The corresponding expression for the background is similar to Eq. ([6]). Adopting the mea- 
sured background flux <p BG (E), which is expected to get much better precision by GLAST 
and PAMELA, and if we know the effective area of the detectors and the identification ef- 
ficiency for the charged particles (of hadronic and electronic background) and photons we 
can get the sensitivity of the detector (l4j. The 'sensitivity' means for some time exposure, 
for example, for one year, the minimal flux the source has in order to have a 5<r detection. 

n 

The sensitivity is not difficult to estimate for EGRET and GLAST, as given in [14J. 
However, a careful simulation of the detector is beyond the present study. We will take 



the sensitivities of EGRET and GLAST directly from the literature 46j, that is, 3 x 10 



and 2 x 10 10 ph cm 2 s 1 respectively. In 46| the sensitivities of EGRET and GLAST 



are for one year of all sky survey with the diffuse gamma background from EGRET as 
2 x 10~ 5 ph cm~ 2 s^ 1 sr^ 1 (lOOMeV/ E) 1A , the typical background at high galactic latitudes. 
Considering background at different latitudes and longer obervation time will certainly 
change the sensitivity. However, the exact values of the detector sensitivities are not very 
important, since the constraints on the SUSY parameter space given in the following can be 
simply rescaled with the sensitivity. From Fig. [2] we can easily understand this: if sensitivity 
is lowered by a factor n the particle factor can be probed is also lowered by the factor n. 

For GLAST we have calculated the similar result to that in Fig. [2] with better angular 
resolution. For Moore profile there is very small difference from Fig. [2j which means most 
annihilation takes place within the very small innermost region at the halo center. For the 
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NFW profile there are difference for the number of the brightest 7-ray sources, which may 
be from subhalos near the Sun. For these nearby sources different angular resolution leads 
to different 7-ray fluxes when taking the NFW profile. 



III. CONSTRAINTS ON THE SUSY PARAMETERS 

From the results in last section we can predict the number of 7-ray sources detectable in 
EGRET or GLAST for any SUSY models. Conversely, if no source is detected the SUSY 
models are constrained. 

We will work in the frame of MSSM, the low energy effective description of the funda- 
mental theory at the electroweak scale. By doing a random scan we give how the parameter 
space is constrained by these detectors. 

However, there are more than one hundred free SUSY breaking parameters even for the Im- 
parity conservative MSSM. A general practice in phenomenological studies is to assume some 
simple relations between the parameters and greatly reduce the number of free parameters. 
Following the assumptions in DarkSUSY [47] we take seven free parameters in calculating 
dark matter production and annihilation, i.e., the higgsino mass parameter /z, the wino mass 
parameter M 2 , the mass of the CP-odd Higgs boson m^, the ratio of the Higgs Vacuum 
expectation values tan/3, the scalar fermion mass parameter mj, the trilinear soft breaking 
parameter A t and A/,. All the sfermions have taken a common soft-breaking mass parameter 
m^; all trilinear parameters are zero except those of the third family; the bino and wino 
have the mass relation, Mi = 5/3 tan 2 6wM 2 , coming from the unification of the gaugino 
mass at the grand unification scale. The simplification of the parameters actually does not 
decrease the generality of our discussion, since the seven parameters are the most relevant 
ones for our purpose. Including other parameters will not change our results much. 

We perform a numerical random scan in the 7-dimensional supersymmetric parameter 
space using the package DarkSUSY [47]] . The ranges of the parameters are as follow- 
ing: 50GeV < \fj,\, M 2 , M A , m f ~ < lOTeV, 1.1 < tan/3 < 61, -3m, < A t , A b < 3m 9 ~, 
sign(/i) = ±1. The parameter space is constrained by the theoretical consistency require- 
ment, such as the correct vacuum breaking pattern, the neutralino being the LSP and so 
on. The accelerator data constrains the parameter further from the spectrum requirement, 



the invisible Z-boson width and the branching ratio of b — > sj 471 ]. 
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FIG. 3: The points in the figure represent models produced randomly in the SUSY parameter 
space. The narrow strip represents the critical values of (av) for these models, as explained in the 
text. In the following figures we use thick curves to represent the strip. 



The SUSY models are divided into two groups: those satisfy the constraint of dark 



matter relic density within 4<r for Q x h 2 



u.ruo_ a013 



48( | and those do not satisfy, i.e., 



Q v h 2 < 0.053. The effect of coannihilation between the fermions is taken into account when 



calculating the relic density numerically. For the second 



group of models we assume the 



neutralino is produced by some nonthermal mechanism 49| to satisfy the observation. 

We then derive the constraints on the SUSY parameter space. Having known the sensi- 
tivities, according to the result in Fig. [2] we can get the detectable number of 7-ray sources 
in EGRET and GLAST for the SUSY model we take in last section. When we scan in the 
SUSY parameter space, we calculate the average detectable number of subhalos at EGRET 
and GLAST for each SUSY model. Then we scale the value of (av ) to a critical value that 
only one subhalo can be detected. The SUSY parameters with larger (av) than the critical 
value of (av ) should be excluded by the experiments if null results are gotten. In Fig. [3] 
we show the critical values of (av ) according to the procedure above, which form a narrow 
strip. We would expect that the constraints should be divergent since the annihilation final 
states should be very different. However, according to Fig. [3]we know that the 7-ray spectra 
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FIG. 4: Constraints of EGRET on the SUSY parameter space if no 7-ray sources are detected at 
the 5a level from dark matter clumps. The left panel gives the constraints assuming a NFW profile 
while the right is for Moore profile. The points in the figure represent models produced randomly 
in the SUSY parameter space. Models above the curves are ruled out. Different curves are given 
adopting different simulation results, as explained in the text. The models which satisfy the relic 
density within 4<r (green points) have smaller (av) than these having £lh 2 < 0.053 (red points). 
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FIG. 5: Same as Fig. |4] except that the constraints are from GLAST. 



from different final states with same m x should be quite similar so that we get convergent 
result. In the following figures we use thick curves to represent these strips. 

In Fig. H]we show the constraints of EGRET on the SUSY parameter space if no uniden- 
tified 7-ray sources are from dark matter clumps. The left panel gives the constraints 
assuming NFW profile while the right is for the Moore profile. We notice for NFW profile 
only the 'subhalo' scenario can constrain the non-thermal SUSY models. For Moore profile 



13 



the 'subhalo' scenario can also constrain a part of the thermal SUSY models. The other 
scenarios have quite weak constraints on the models by EGRET. For light neutralinos the 
constraints on (av) are 10 _26 cm _3 s _1 and 10~ 27 cmr 3 s _1 for the NFW and Moore profiles 
respectively in the 'subhalo' scenario. For other models the constraints are ~ I0~ 24 cm~ 3 s~ 1 
and ~ 10 _25 cm _3 s _1 for NFW and Moore profiles respectively. We notice that the models 
with Qh 2 < 0.053 have a greater (av) than those thermal models and therefore produce 
larger 7-ray fluxes. These models are easier to be ruled out. 

In Fig. [5] we show the similar constraints on SUSY by GLAST, which can give a much 
severer constraints on the parameter space than EGRET gives. For the Moore profile all 
scenarios can put constraint on the SUSY models. In this case a large fraction of the 
nonthermal parameter space will be ruled out. The constraints on (av ) for light neutralinos 
are now 10 _28 cm~ 3 s _1 and 10 _29 cm _3 s _1 for the NFW and Moore profiles respectively in 
the 'subhalo' scenario. For other models the constraints on (av) reach ~ 10~ 26 cm~~ 3 s~ 1 and 
~ 10~ 27 cm~ 3 s _1 for NFW and Moore profiles respectively. 



IV. DISCUSSIONS AND CONCLUSION 



Since the rate of dark matter annihilation is proportional to the DM density square, the 
Galactic center had been considered as the most promising site to search for the annihila- 
tion signals. The possibility of detecting dark matter annihilation from the GC has been 
extensively studied in literature However, the GC is a very complex environment. The 
dark matter density profile near the GC is complicated due to the existence of baryonic 
matter and leads to difficulties in making theoretical calculations. For example, the SMBH 

can either steepen or flatten the slope of the DM profile at the innermost center of the 

I I 

halo depending on the evolution of the black hole [50J- Furthermore, the baryonic processes 
associated with the central supermassive black hole (SMBH) and the supernova remnant 
Sgr A* [51] provide a strong 7-ray background, which has been detected by H 
the signals of dark matter annihilation and make the detection very difficult 



SS m^to 



5l|. In [51| 



it is found that in case of NFW profile GLAST can probe (av ) between 10 



-26 



cm 



1 and 



10 cm s for light neutralinos, which is similar to the sensitivity by observing subhalos 
as shown in Fig. However, considering the HESS detected 7-ray background only models 
with (av) > 1§~ 27 cm~ 3 s~ x can be probed 5l| from the GC by GLAST. Therefore it becomes 
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less sensitive to probe the GC than detecting subhalos now. For the case of Moore profile, 
sensitivity from GC observation is improved by two orders of magnitude, while only one 
order of magnitude improvement from subhalo observation. This means similar sensitivities 
will be achieved by observing the GC and subhalos. 

On the contrary, subhalos provide a clean environment to search for the annihilation 



signals. Especially recent simulation shows that t 



45 



re DM profiles may not be universal. 
531 ] . Reed et al. gave the cusp index 



Smaller subhalos may have steeper central cusp 
7 = 1.4 - 0.081og(M/M*) for halos of 0.01M* to 1000M* with a large scatter. In this case, 
if taking the GC the NFW profile and the subhalos the Moore profile, the 7-ray fluxes from 
the subhalos may even be greater than that from the GC. 

We expect these sources can be detected by the satellite based experiments, such as 
EGRET and GLAST. Once such sources are detected we can learn a lot about the nature of 
dark matter particles by studying its luminosity and spectrum of the annihilation. However, 
study shows that most of the EGRET unidentified 7-ray sources should not be of dark 

nn . 

matter origin |6|, 123[ . In this work we study how the EGRET and GLAST can constrain the 
SUSY models if none of the subhalo 7-ray sources are detected in the two experiments. 

We first realize one hundred MW-like halos with subhalos whose distribution is given 
according to the N-body simulation results. In each realization we calculate the 7-ray 
flux from the subhalos by fixing the particle factor. Then we give the average cumulative 
number of the subhalo 7-ray sources as function of their flux intensities. Once the sensitivity 
of detectors, such EGRET and GLAST, are known, we know the detectable number of 
this kind of 7-ray sources. By requiring the detectable number smaller than 1 at EGRET 
and GLAST we put a constraint on the SUSY parameter space. Our result shows that 
the EGRET has already given a moderate constraint on the SUSY parameter space if we 
assume none of the unidentified 7-ray sources are from subhalos. The GLAST can greatly 
enhance the constraints. However, a large uncertainty comes from the simulation, especially 
the property of subhalos in a dense environment. Convergence of the subhalos property in 
the future high resolution simulation will lead to more precise constraint on the nature of 
dark matter particles. 
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